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Abstract 

We prove a meromorphic integrability condition at order 2 near a homothetic 
orbit for a meromorphic homogeneous potential of degree —1, which extend 
the Morales Ramis conditions of order 1. Conversely, we prove that if this 
criterion is satisfied, then the Galois group of second variational equation is 
abelian and we compute explicitly the Galois group and the Picard-Vessiot 
extension. 
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1. Introduction 

We want to study the Hamiltonian systems of the following form 

H^T + V 

2 

with 7" = X] -f and V a homogeneous function of degree — 1 in gi, . . . , gf„ and 
n e N*. The corresponding differential equation system is written 

dpi dqi 

For the following and our non-integrability proofs, we will only consider mero- 
morphic potentials V . 
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Definition 1. We call c G C" Darboux point a solution of the equation 

V'{c) = ac 

with a G C called multiplier. Because of homogeneity, we can always choose 
a = 0, — 1, and we will say that c is non degenerated if a 7^ 0. 

Definition 2. A homothetic orbit associated to a Darboux point c with 
multiplier a is given by 

qi{t) = Ci(j)(t) pi{t) = Ci4>(t) i = l...n 

with (j) satisfying the following differential equation 

1 , ; n2 a ^ 

with E corresponding the energy orbit. 

In the following, we will call "norm" and scalar product the expressions 

n n 
|f P=^^fj^ < V,W >= VjWi 

i=l i=l 

even for complex v,w. We will say moreover that a matrix is orthonormal 
complex if its columns Xi, . . . , X„ are such that 

n n 

<X„X, >=^(X,),(X,)fc = Vz,j \X,\^=Y,{Xi)l = l Vz 
fc=i fc=i 

Theorem 1. (Morales, Ramis, Yoshida 11^,1^,13^,1^) Let V be a meromor- 
phic homogeneous potential of degree —1 and c a non degenerate Darboux 
point. IfV is meromorphically integrable, then the identity component of the 
Galois group of the variational equation near the homothetic orbit associated 
to c is abelian at any order. Moreover, if we fix the multiplier of the Darboux 
point c to —1, then the identity component of the Galois group of the first 
order variational equation is abelian if and only if 

Sp(VV(c)) C |^(A; - l){k + 2), ke n\ 



2 



The main theorems of this article are the following 

Theorem 2. Let V be a meromorphic homogeneous potential of degree —l,c 
a Darboux point ofV with multiplier —1. We will suppose that V'^V{c), 
the Hessian of V in c, is diagonalizable. We pose Aj i — 1 . . .n 
its eigenvalues and Xi, . . . , its eigenvectors. If V is meromorphically 
integrable, then 

J _ (Pt-iKp»+2) y^^^^ pj e N (first order integrability condition) 
ii \/i,j,k = l...n, Ap^.p^^p^ = ^ D'^V{c).{Xi,Xj,Xk) = where A is a 3 
index table with values in {0, 1} invariant by permutation and given by 

• For i,j,k G N*, Aij^k — ^ if (^nd only if one of the following condi- 
tions are satisfied 

-i+j + k<-3 
i + j + k mod 2 = 1 

f i-j + k<-3 ( i + j-k<-3 

\ i + j + k mod 2 = 1 \^ i + j + k mod 2 = 1 

• For i — 0, j,k e N*, Aoj,fe — 1 if and only if \ j — k |> 2 

• For i = j = 0, Aij^k = 1- 

Theorem 3. Let V be a meromorphic homogeneous potential of degree —1, c 
a Darboux point ofV with multiplier —1. We will suppose that V^l^(c), 
the Hessian of V in c, is diagonalizable. If V is meromorphically 

integrable, then the Galois group of the 2-th order variational equation is 
always isomorphic to C and the Picard- Vessiot extension field is 

c(^0,0,lnQ + 0(^l + -^0 

except if one (or both) of the two following conditions are satisfied 

• D^{V){c){v,v,v) ^ with VV(c)v ^-v, v^O 

• D^{V){c){v,v,w) ^ with VV(c)^; = -v, VV(c)w = 0, v.w^O 



i+j-k>2 
i-j+k>2 
-i + j + k>2 
i + j + k mod 2 — 



or 
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and in this case the Galois group is and the Picard-Vessiot field 



2 V V2' 

The first order condition is already known, computed by Yoshida P based 
on classification of hypergeometric functions by Kimura [5] . It has been used 
many times by [6], [7], in particular in the n body problem in the case of 
homogeneity degree —1 in [5]. The Morales Ramis theorem holds for varia- 
tional equations at any order, and so here we want to study completely the 
second order, and give an integrability characterization at order 2. The inte- 
grability constraint, when the system of second order variational equations 
is well written, can be found by computing particular monodromy commuta- 
tor, like done in [Hj. But the true difficulty is not here, it is that we need to 
study the monodromy for an infinite number of eigenvalues, for all possibles 
eigenvalues. It comes down to the study of a particular 3 index sequence, and 
then to find its zero and non zero entries. Of course it can be easily checked 
one by one, but this is not enough. By chance, this 3 index sequence possess 
an explicit expression, but not easy to find (and to prove). The property 
behind it is that the monodromy commutator is D-finite with respect to the 
eigenvalue parameters, and so it satisfies a 3 index linear recurrence with 
polynomial coefficients. A closed form solution can then be guessed by the 
gfun package, and then its validity checked. The non nullity can then easily 
be studied because this closed form expression is a hypergeometric sequence. 

2. Variational equations of order 2 

2.1. Reduction of second order variational equations 

For the following, 0(t) will be a solution of the equation = ^ + 1. 
The second order variational equation, if the Hessian matrix is diagonal, can 
be written 

1 1 / yityT.Yit) \ 

X = ——DX + — ■■■ Y = ——DY (1) 

m 2 0(t)4^^^^),^^^^^J 0(t)3 

with D diagonal, Tj G M„(C). The matrix D is the Hessian matrix of 
V on the Darboux point c, and the matrices Tj, i = l...n are defined 
by Tjjfe = D^(y){c).{qi,qj,qk). A more detailled construction of second 
variational equation and higher orders can be found in [9]. 
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Theorem 4. We consider the differential equation ([T]). Suppose that the 
equation in Y has a virtually abelian Galois group. Let us note K the Picard 
Vessiot extention of ([T]), generated by all solutions of Q. The Galois group 
of ([!]) is virtually abelian if and only if all the equations 



have virtually abelian Galois groups. Moreover, the Galois group depends 
only on the nullity or non nullity ofTi j^. 



Proof. Suppose that the Galois group of Q is virtuaUy abelian. We write 

Clflit)+C2f2it) 
C2n~lf2n-l{t) + C2nf2n(t) 



Y{t) 



using the fact that D is diagonal. Let us fix all the Cj = except C2j-i, C2j 
for some fixed j. We come down to the equation 

2 11/ -^I'i'i^'^^)^ 



X{t) = -—rDX{t) 



<P{tf ' ' 2<p{ty 



2 



The function Yj{t) is a fixed function, so it is a non homogeneous linear 
differential equation. By hypothesis, the Galois group of this equation is 
virtually abelian. Using the method of variation of the constant, we find that 
the Galois group depend only on the nullity or non nullity of this equation 

We can then substract to T all its terms of the form Tj and still the 
Galois group of ([T]) will be virtually abelian. Now we do the same procedure 
for the terms in Yk{t)Yj{t). We fix all the C = except for C2j-i,C2j and 
C2A;-i, C2k- As we have Tijj = and Ti^k,k = we get the equation 



X = -^DX+ ^ 



<p{tf 2<p{ty 



T,,^kY,{t)Yk{t) 



The functions Y^(t), Yfc(t) are fixed, so this is a non homogeneous linear dif- 
ferential equation. With hypothesis, the Galois group of this equation is vir- 
tually abelian. Using the method of variation of the constant, we find that 
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the Galois group depend only on the nullity or non nullity of this equation 
Tijj. Conversely, the equations ([2]) are non homogeneous linear differential 
equations, so we can sum the solutions of the equations 

Xi = ^^^Xj + TijkYj(t)Yk{t) 

to produce all the solutions of equation ([T]). 

The following lemma will have a primary importance in computation of 
monodromy. In fact, it will be necessary only to compute some sort of residue 

Theorem 5. We consider F G Clzi) [Z2] and 

fit) = F ^t, arctanh ^- 
We consider the differential field and the Galois group 

K = C arctanh ( ^) , / M) G = a{K,C{t)) 



If G is ahelian, then 



■^Res F arctanh + = Va G C (3) 



Proof. First we recall that if the Galois group G is abelian, then so is the 
monodromy group, because the monodromy group is always included inside 
the Galois group. We consider two paths, the eight path ai around the 
singularities —1, 1 and the path around (T2. At infinity, F (t, arctanh (^) + a) 
has the following series expansion 

J F ^t, arctanh ^-^ + dt = a„(a)t" + r{a)ln t 

because the function arctanh is smooth at infinity. We consider now the 
commutator 
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X singularities sigma_2 sigraa l 



Figure 1: Paths corresponding to monodromy elements ai,a2, and the Riemann surface 
associated to the Qj. The difference between two sheeves is inejPj. We see that (72, 
corresponding to monodromy around infinity, acts triviaUy on Qj . 



We have that (xf " (/) = F (t, arctanh (|) + /3) and a2(/n t) = lnt + 2in. We 
conclude that 

^(/) = / + r(/3)-r(0) 

This r(a) correspond to the residue of F (t, arctanh (^) + a) at infinity. If 
the monodromy is abelian, then a should act trivially on /. This is the case 
if and only if r(/3) — r(0) V/3 G Z. The function r is polynomial in /3, then 
r{l3) - r(0), V/3 G C. This gives the formula Q 



2. 2. Computation of terms of order 2 

Now we prove that we can always come down to study the equation ^ 
to prove Theorems [2l [3] 



Theorem 6. Let A G M„(C) be a symmetric matrix, which mean Aij = 
Aji. We suppose that A is diagonalizable. Then A is diagonalizable in an 
orthonormal basis, it means that it exists Xi, . . . such that 

n n 

<X„X,>=Y,{X^)k{X,)k = I X, |2= = 1 Vz 

k=l k=l 



7 



This small result is in fact very important for reduction of problems deal- 
ing with homogeneous potentials. This is because the integrability status of 
a potential V is not changed after a rotation (even a complex one) nor dilata- 
tion. Such property allows a great simplification in classification like in [6], 
[7]. The simplification at order 2 is even more important, as a such simple 
criterion of Theorem |2] is not possible if the Hessian matrix is not diagonal. 
As we will see, the only hypothesis of diagonalizability is very weak, because 
in particular the system is rarely integrable at order 1 if the Hessian matrix 
is not diagonalizable (the conditions on the spectrum are still necessary, but 
not sufficient). 

Proof. Let A G M„(C) be a symmetric diagonalizable matrix, Xi, . . . ,X„ 
and its eigenvectors. Suppose first that | VA;. We have then 

< Xk, Av >= ^(Xfc)i ^ aijVj = X] = 

i j j i 

< A^Xk, V >=< AXk, v>= Xk< Xk, V > 

Then the orthogonal subspace to X^ is stable by A. Moreover, is not 
in Ek then E^ is supplementary to CX^, then E^ containts all the other 
eigenvectors. Then 

n 

<X„X,>=Y,iX^)k{X,)k = Vz,i 

k=l 

Now suppose it exists X^ such that | Xj. P= 0. Then X^ is orthogonal to all 
eigenvectors with non zero norm (using the proof before). So if it is the only 
eigenvector with zero norm, then Ek containts all the other eigenvectors, but 
also containts X^ because < X^^X^ >= 0. But E^. is of codimension 1, then 
it cannot containt all the eigenvectors (because form a basis). 

Then it exists another vector, let say Xi, with zero norm, and moreover 
< Xk,Xi >7^ because otherwise the previous reasoning still holds. 
Let us now write the matrix P G M„(C) formed by the eigenvectors Xi, . . . , 
Xn- We have then 

A = PDP-^ = = [P-^DP^ 
where D is a diagonal matrix. Then 

P^PD = DP^P 



8 



Then P^P commute with D. We know that [P^P]^^^ ^ 0, then Dk,k = Di,i 
(because P^p should commute with D). Then Xk and Xi have the same 
eigenvalue. So, the fact that these vectors have zero norm only correspond 
to a bad "choice" of eigenvector. It suffice to take + Xi and X^ — Xi 
instead of X^, X/. 

Now that I Xfc Vfc, we just need to multiply X^ with some constant 
such that I Xk P= 1 VA;. This end the proof of the theorem. 

Thanks to this theorem, we know that if the Hessian matrix is diagonal- 
izable, we can always make an orthonormal variable change in the potential 
V (this does not affect at all integrability) such that the Hessian matrix 
becomes diagonal, and thus apply Theorem |4j Let us check now that the 
condition of Theorem [2] is well defined and equivalent to this. 

Theorem 7. Let V he a meromorphic homogeneous potential of degree — 1 
and c G C" a Darhoux point with multiplier —1. We suppose V^l^(c) diag- 
onalizable and we note its eigenvectors Then the integrability 

constraints of Theorem^ do not depend on the choice of Xi, . . . ,X„. More- 
over, if we make an orthonormal choice for the Xi, . . . , X„, then coefficients 
Tij^k in ([T]) will correspond to D^(y).{Xi^Xj^Xk). 

Proof. First, if we make the orthonormal choice (which is always possible 
thanks to Theorem |6]), we note P the associated orthonormal matrix. We 
have then that the potential W{q) = V{Pq) has a Darboux point in P~^c 
and the corresponding Hessian matrix is diagonal. Moreover, we have that 

D\W).{q„q„qk) = D\V).{Pq,, Pq,, Pqu) = D^(\/).(X„ X„ X^) 

Now, we verify that the criterion is well defined. We first consider the case 
where all eigenvalues of V^K(c) are distinct. Then, up to multiplication by a 
non zero constant, there is a unique choice of eigenvectors Xi, . . . , X„. So this 
does not change nullity or non nullity of {V) . (Xj , Xj , X^ ) . Now, if there are 
multiple eigenvalues. There are an infinite number of choices for eigenvectors 
X. We fix one. Suppose that Xi,X2 have the same eigenvalue. Since the 
nullity condition is associated only to the corresponding eigenvalues, if there 
is a nullity condition for some third order derivative involving Xi, it will be 
the same for X2. Suppose there is a condition 

D\V).{X^,X,,Xk) = D\V).{X2,X,,Xk) = 
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Then, for the vector aXi+pX2, we will also have D^(y).{aXi + l3X2, Xj,Xk) 
= by expanding it. Remark that even if j = 1, A; = 1, it will still work. All 
basis changes can be written as successive such linear combinations. So the 
constraints of Theorem |2] do not depend on the choice of Xi, . . . , X„. 

This theorem is in principle not very important, because we can reduce 
the potential directly to have a diagonal Hessian matrix, but it is useful 
in pratice. This is because it is not so easy to compute the orthonormal 
basis change as just to compute the spectrum. In particular when there 
are parameters, the normalization of eigenvectors becomes uneasy because 
this produces singularities. We just say here that we can avoid that, and 
do not care about normalization (for some specific values of the parameters, 
Xi,. . . ,Xn will no more be a basis, but the criterion is in this case still a 
necessary one). 

3. Non integrability of second order variational equations 

3.1. A first approach 

Integrability condition at order 2 for equation ([T]) 

Theorem 8. We consider an equation of the form 

with D diagonal, Tj G Mn(C). The equation ([T]) has a virtually abelian Galois 
group if and only if 

• Di^i = with Pi (z'R ( integrability condition of order 1 ) 

• Wi,j, k = l...n, Ap.^p.^p^ = ^ Tij^k = where A is a three index table 
with values in {0, 1}, invariant by permutation and whose first values 
are given by 
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One direct application of this theorem is the study of problems in celestial 
mechanics, like in [8], [in]. It is often unnecessary to know the table A for 
arbitrary high eigenvalues, except in some cases like the open problem at the 
end of [7j, which is because of this much more difficult. 

Proof. Using Theorem |4| the study of the Galois group of ([T]) comes down 
to the study of Let us look first at the expressions of functions Y. We 
can fix E = 1 because we can fix E to any value using homogeneity (except 
for E = 0). The variable change (f){t) — > t gives the equation 

We then make the variable change 



1 + t 



t 

The equations satisfied by Yi{t) becomes 

1/9 X - 1) + 2) 
-{f-l)y + 2ty-^- ^-j '-y = Q 

A basis of solutions is given by (Pj, Qj) where Pi are polynomials and the 
functions Qi can be written 

The functions Q are multivalued except for i = 0, which is particular, because 
the Galois group is Id instead of C and then all solutions are algebraic. We 
get for X{t) the following solution 

X{t) = CiP.it) + C2Q^{t) + 

j Yj{t)Yk{t)P,{t){t^ - lfdtQ,{t) - J Y,{t)Yk{t)Q,{t){t^ - lfdtPi{t) 
So we need to study the monodromy of 

j Y,{t)Yk{m{t){e-lfdtQ,{t)- j Y,{t)Yk{t)Q,{t){t^-lfdtPi{t) 
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Let us try to apply Theorem |5} This theorem do not apply directly because 
there could be compensations between the two integrals. But we can rewrite 
it 

QAt)Qk{t)m){t'-ifdtp,{t) J ^^^—L^dt- 




Qjmkit)m)dt^^^-^^j^^it' - ifdtpm 



Then 



/ 



Qj{t)Qkit)P,{t)dt 



IS m the Picard Vessiot field of Q (b ecause Pi is a polynomial and that the 
Picard Vessiot field is stable by derivation). We also have that Qi is in the 
Picard Vessiot field, and then by substraction, 

Qjit)Qkit)Qi{t){t''-lfdt 

IS m the Picard Vessiot field of Q. We can now apply Theorem |5] to this 
integral, and so it is only needed to study the residue 

S = Res if - lf{Q^{t) + eiaPi){Qj{t) + ejaPj){Qk{t) + ekaPk)dt 

t=oo 

The polynomials Pi can be generated by the formula 

(which gives a normalization for the dominant coefficient of Pi that we will 
choose for now) and the functions Qi can be written 



Qi{t) = eiPiit) arctanh 



l\ , W,{t) 



t t'-l 



with Wi polynomials, a real sequence. 

The integrability at order 2 only require that the identity component of 
the Galois group be abelian. This does not say a priori anything about the 
whole Galois group, except if it is connected. This is the case here. The 
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Galois group at order 1 is C (or Id). Then at order 2, the Picard Vessiot 
field will be of the form 

K = c(^, arctanh , J f{t)d?j with f{t) e C arctanh ^ 

We just add some integral in K. Then the Galois group (T{K,C{t)) is still 
connected. 

So, we write in the table Aij^k, 1 is the constaint of Theorem |5] is satisfied, 
otherwise. This criterion is a priori only a necessary criterion, not sufficient. 
So we then try to compute the integral using integration by parts 



QAt)Qk{t)mW - iydtQ,{t) - J Qj{t)Qk{t)Qi{t){t' - lydtPiit) 

Using the expressions of functions Q, we need to integrate functions in 
C [t, arctanh {]:)]■ We make successive integration by part, by deriving the 
term in arctanh of highest degree, and finally we arrive on an integral 
of a function in C(t). This procedure could fail, but it works every time 
for Aij^k = 1 (for Aij^k = 0, the procedure fails because terms in ln(t^ — 1) 
appear) . Using Theorem |4| we know that it is necessary and sufficient that 
all equations 

have a virtually abelian Galois group for all non zero Tj ^ for the virtual 
abelianity of the Galois group of ([!]). 

Theorem 9. The table A of Theorem^has the following values 

• For k G N*, Ajj,fc = 1 ^/ (I'fT'd only if one of the following conditions 
are satisfied 

i+j-k>2 

i-j + k>2 ( -i + j + k < -3 



or ^ 

-i + j + k>2 \ i+3 + k mod 2 = 1 

i + j + k mod 2 = 

f i-j + k<-3 ( i + j-k<-3 

\ i + j + k mod 2 = 1 \ i + j + k mod 2 = 1 

For i = 0, j,k &W , Ao,j,fc = 1 ^/ (iiT'd only if \ j — k \>2 
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• For i = j = 0, Aij^k = 1- 

Moreover, the table A is invariant by permutation of the index i,j, k. 

This table is the direct generahzation of the integrabihty table of [T] at 
order 2 for degree —1. A similar process could be done for other homogeneity 
degrees, but in fact the degree —1 is much more simple for three reasons 

• There is only one family in the Morales Ramis table for degree —1, and 
generically, there are two (and the complexity increase with the power 
three of the number of families) . 

• Some homogeneity degrees have very particular families, associated the 
groups A4, 5*4,^5. This produce very very complicated computations. 

• By studying only one homogeneity degree, we have one parameter less 
in the D-finite computations. This is important because computational 
cost usually increase exponentially with the number of parameters (at 
least). 

3.2. Study of the solutions 

Proof. Using the last theorem, we already know that we just have to study 
the residue 

Res - l)\Q,{t) + e,aP,)(Qfe(t) + tkaPk){Qk{t) + e^aPfc) 

and a necessary condition for integrabihty is that this residue should be 
independent of a. We will call the fact that the coefficient in should be 
zero "the constraint in a^" and respectively "the constraint in a" for the 
term in a. We see also that a priori, the residue is a polynomial of degree 3 
in a but we have that 

coeff (Res {Qi{t) + e,aP,){Qj{t) + e,aP,)(Qfc(t) + ekaPk){t' - l)',a=^) = 

\t=oo J 

Res P,P,Pfee,e,efc(t2 - 1)^ = 

t=OD 

because it is a polynomial in t. So we have only two constraints for inte- 
grabihty. Let us begin by checking that all functions Qi are multivalued for 
i 7^ 0. We only need to prove 

Lemma 10. We have ti 7^ Vi G N* 
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Proof. Looking at first values of ej, we already can guess the expression of 

^ 4-H(i + 1) 

We now need to prove it. The sequence can be computed thanks to the 
formula 

f 1 

:dt 



I it' - ly 



c 

with C a circle around —1,1 in the direct way (because ej is the term in front 
of arctanh which grows by 1 along C). Using the symmetry t — )■ — t, we 
only need to compute the residue in 1 for example. We have then 



d 



\t=\ 



dt \(t + l)^Pi{t)\ 

knowing that 1 is never a root of Pj. So we just need to compute the sequences 

p,(i), (J-^m)^ \t=i 

with the recurrence formula Pi 

(4n2 + 12n' + 8n)P„ + {-Atn'^ - Utn - 12t)P„+i + (n + 3)P„+2 
We have then 

P,(l) = + 1)! \t=i= + 3){t + 1)! 

and we get e^. 

For the following, we will use the system 

{{4n^ + Vlr? + 8n)/„(i) - (4in' + 14tn + 12t)fn+i{t) + + 3)/n+2(0, 

- l)m + Atm -(n-l)(n + 2)Ut)} 

which vanish for fn{t) — Pn and fn{t) = e~^Qn- The system 

{-4tut) + it' - -ut) + 

vanish for = (t^ — 1)^. We will use these systems for the package Mgfun 
to compute recurrences for our residues. 
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3.3. General case 

Proof. First part. We prove that the Galois group is not virtually abelian 
if the variational equation contains a term corresponding to index such that 
■^i,j,k = 0. We will begin with the non zero index case (this is because the 
index is very special, in particular the function Qq is not multivalued). We 
now need to compute the residues of Theorem |5] for all index and prove they 
are non zero for = 0. Knowing that 7^ i > 1, it comes down to the 
study of the sequence 



Si,,,k = Res (er'Qiit) + aPi){e-^Q,{t) + aPj){e^^Qk{t) + aPk){t' - 1)' 

t=oo 



We have moreover that the system Sys vanish for P„ and e'^Qn, and so is 
also vanishing for e~^Qn + ct-Pn- Thanks to that, we will be able to find a 
recurrence on Si^^k-, and moreover it will not depend on a. 

Lemma 11. The sequence Sij^k satisfy the following recurrence relation 



-{l + i){i+j + k + 2)Sij^k + 4^(2 - l)(i - 2)(i - 3 - J - k)Si.2,j,k+ 
4t{2t - - l)^,„i,,_i,, + At{2t - l)k{k - l)S,_i,,,fc_i = 



We have also that Sij^k is invariant by permutation of the index, and this 
recurrence relation is not, and so this produce other recurrence relations. 

This recurrence relation can be proved automatically using the Mgfun 
package for Maple, the holonomic package for Mathematica, or even at hand 
using integration by parts and a formula between the derivative of Pn and 
Pn,Pn-i- Now we will in fact prove much more than necessary for proving 
Theorem [9} we will find closed form solutions for the residue we need to 
compute. 

Lemma 12. For k > j > i > 1, we pose 



with a = —i + j + k, b = i — j + k, c = i +j — k,d = i + j + k. The expression 
of Sij^k is given by 



fihj, k) 



2^.!j!A:!r(i(rf + i))r(f)r(|)r(f) 



r (i(a + 3)) r {\{b + 3)) r (i(c + 3)) r {\{d + 4)) 




Si,j,k = eval{Sij^k, a = 0) + 
i + e, j + e, k + e)a^ if i + j + k mod 2 = 1 
f{i + e,j + e,k + e)a if i + j + k mod 2 = 
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The limit can be easily computed for all i,j,k G N* but there are no 
closed form expression for the limit valid for all i,j, k. The limit depend in 
fact of the order of i,j,k. We choose in fact such a complicated formula 
because of its generality. It holds in all cases and thus allows to speed up 
the proof, avoiding to do 3 times the same thing, and show effectively the 
symmetry between the index. With these formulas, it will be easy to prove 
Theorem [9] because the cases daSij^k = correspond to singular values of the 
r functions in the denominator. 

Proof. 3.3.1. Case i + j + k mod 2 = 1 

We begin by looking at / for i + j + k mod 2 = 1. This is the easy 
case, because when we replace e by in /(i + e, j + e, + e), the expression 
f{i,j,k) i,j,k e N* is still meaningfuU if we assume r(— n) = 00, n e 
N. Indeed, there can be at most one term of this kind and always in the 
denominator. The corresponding value of Sij^k will be 0. We then check 
than this expression formally satisfy recurrence (|4]). First of all, we remark 
that when we select the coefficient in Sij^k, we find 

coeff(5'jj_fc, a^) = SRes PjPj arctanh ( - ) (t^ — 1)^ 

t=oo V ^ / 

Then we make a series expansion of arctanh (|) for t = 00, and we also get 
that (noting that PiPjP^ is an even polynomial) 

p=o 

This relation is interesting because we have an orthogonality property on the 
polynomials Pj. We now need to study boundary cases. Using the recurrence, 
we come down to i = 1,2. The recurrence relation Q get simpler for i = 2 
and so we can express 5*2 j,a: in function of Sij^k- Using orthogonality, we get 
the formulas 

1 

j PjPkit'' - Ifdt = J - A; ^ -2, 0, 2 

-1 

which give Sij^k = 0, j — k ^ —2, 0, 2. We get for f{l,j, k) = 

48(_l)|0-fc) sin (ivr (4 + j - k)) 2^+^! {k + l)T (j + 1) 
~{2+j-k) {j -k){j-k-2)7r {3+J + k) (j + 1 + k) (j + k- 1) 
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Because of the sinus term, this formula vanish for j — k —2, 0, 2. One just 
need to check it is right in these three cases. We have 

^9.9 , 3 4'=r(A;-l/2)r(fc + l)2 
P^(t^ - ifdt = ^ \ 

-1 

' 2 T{k + 7/2) 

using Mgfun, and we check these are the same as /. 



-1 



3.3.2. Case i+j + k mod 2 = 

We now look at the function / for i+j + k mod 2 = 0. This time, if we 
replace formally e = in p^/(^ + e, j + e, k + e)a, we find a quotient of F 
functions and in the numerator at most a term of the form r(— n), n G N. 
We can still regularize the formula using the relation r(n + 1) = nr(n). We 
get in particular that the limit 

limr^/(2 + e,j + e,fc + e) 

is always finite. If there are no term of the form r(— n), n G N in the 
numerator of /, then the limit is zero. Using invariance by permutation, we 
can suppose that k > j > i, and so the only possible infinite term in the 
numerator is the term in i + j — k. We get then a zero limit for k < i + j, 
and for k > i + j, we can regularize the formula. 

We then check that the formula satisfy the recurrence. Now let us look 
at the boundary cases. With the recurrence, we can reduce i and then Sij^k 
is completely determined by Sij^k, i = 1,2. Looking at relation Q for i = 2, 
it simplifies and we can express S2j,k in function of Sij^k- Using relation Q 
after permutating the index, we get the recurrence 

-(1 + i){i+j + 3)5i,,j + 2){t-j - 4)5i,,_2,i+ 

4i{2i - - l)Si,i-i,j-i = 

which allows us to reduce j to 1. So we only need to compute 5*1, i^a: using 
Mgfun 

lQ''T{k + 3/2)T{k - l/2)r(A; + 2)T{k + 1) 



2+2k 



T {k + 4:)T{k + 5/2) ^/^ 
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Our expression correspond to this one. 



To conclude the proof of table A for non zero index, we now look to the 
expression of /. Looking at the formula of S for i + j + k mod 2 = 1, we see 
that it vanishes exactly when one of the three quantities 

-i+j+k+3 -i+j+k+3 -i+j+k+3 

is non positive. This exactly correspond to the formulas of table A for i + 
j + k mod 2 = 1. In the case i + j + k mod 2 = 0, the quantity 

lim ^/(i + e,j + e,k + e) 

vanish if and only if all the numbers + k, —i+j + k,—i+j + k are 

positiv, which is equivalent using the parity condition to 

-i+j + k>2 -i+j + k>2 -i+j + k>2 

3.3.3. Case of a zero index 

We now look to the case with at least one zero index. We can write 

t 1 /1\ W(t) 

^oW = T^-r Qo(t) = 7^-r Q^(^) = e^P^(^)arctanh' ^ ■ '^'^ 



(the notation Pq and Qq is arbitrary here because they are both rational). 
Wc begin by the case where exactly one index is zero. We need to compute 
the residues 

Res (er^Qiit) + aPi){e-^Qj{t) + aPj){e - 1) 

t=oo ■' 

Res {er^Qiit) + aPi){eT^Qj{t) + aPj)t{t^ - 1) 

t=oo ■' 

These are polynomials of degree at most 2 in a but the coefficient in 
is always zero because we take the residue of a polynomial at infinity. So 
one just need to compute the residue in a. We expand and suppress the 
polynomial terms and we get the formula 

Res (ciej^ + ^j^J^) arctanh ( - j PiPj{t^ — 1) 
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Knowing that e^e^ ^ + CjCj ^ does not vanish, we can suppress it and we only 
need to compute the following sequences 

Slj = Res arctanh ^-^ PiPi+j{t^ — 1) 



5"? • = Res arctanh I - J PiPi^jt{t'^ — 1) 



, = Res arctanh I - 

t=oo \ t 

One just need to prove that either 5*^?^^. or Sfj is not zero for i E N* , j = 0,1 
(the condition on the index of table A such that A = correspond here to 
—2<i—j<2 and we use symmetry of the index). Using the parity on 
if: of the polynomials Pi, we find that only S}q, Sfi can be non zero. These 
sequences can be easily computed for finding a recurrence with Mgfun and 
then a closed form 

3 4T(^ + 1)^ ^2^ = -4- 4T(2 + 1)2 



{i + l){2i + l)i (2i + l)(2i + 3) 

These expressions do not vanish. 

3.3.4- Integrahility in the cases where Aij^k = 1 

Second part: We now prove that if all the non zero terms of second order 
variational equation correspond only to cases such that Aij^^ = 1, then the 
Galois group is abelian. We use the following lemma 

Lemma 13. We consider F{t) = 

iV „ , . , ^lV „ , . , ^l 



H^{t) arctanh |^-J + H2{t) arctanh |^-J + Hi{t) arctanh |^-J + Hi{t) 

with Hi, H2, H3, Hi E C[t]. If the conditions of Theorem^ are satisfied, then 

• If Res F{t) = 0, then J F dt e C [t, arctanh {])] 

• // Res F{t) ^ 0, then J F dteC[t, arctanh (|) , In {t^ - 1)] 

Proof. We proceed using integration by parts. We derive the term in 
arctanh (|) . Posing 

t 

J{t) = [ H,{s)ds 



21 



we have that J(l) = using the condition in o? of Theorem |5] and making a 
series expansion at infinity of arctanh (^) . Then {t^ — V) divide the polynomial 
J. After integration by parts, we get a term of the form 



— 3i?(t) arctanh 



with R a polynomial. Let us try another integration by parts. We get the 
term 

-^TZl J -^^^ + H2{t)dt arctanh (^^ 

We want that this term can be written Z(t) arctanh (|) with Z a polynomial 
(with a good choice of integration constant). We only need that 

1 

-^ + H,{t)dt = 

-1 

Let us look now at the coefficient in a of the residue ^ . We know it is equal 
to zero. 

coeff ( Res F [ t, arctanh (-)+«),« 



(Res F It 

\t=oo y 



Res 3-^3 (t) arctanh - + 2H2(t) arctanh 

t=oo \t J \t 

which gives using an integration by part (we can see the residue as an inte- 
gration along a small circle around infinity) 

Res — ^'^^^^ arctanh ( - ) + 2iJ2(^) arctanh f -] dt 



t=oo t2_X 

Using the Taylor expansion of arctanh (|) at infinity, we get 



1 



-1 

This is exactly our condition (|3]). So the last remaining integral to compute 
is of the form 



Z{t) arctanh [j]dteC 



t, arctanh ^-^ , In (t^ — l) 
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which can be proved using an integration by part. Now let us look closer to 
the possible terms in In (t^ — 1). Let suppose that is exists a term In (t^ — 1) 
in J F dt. We will have 

j Fdt = Z^{t) arctanh +■■■ + Zo{t) + r In (t^ - l) 

with Z^, . . . ,Zq polynomials and r is a constant because In (t^ — 1) do not 
appear in F. A function in C [t, arctanh (^)] is meromorphic near infinity. 
Derivating this expression will give 

with g a meromorphic function on a neighborhood of infinity. Then 

Res F{t) = r 

t=oo 

So, if this residue is zero, there will be no In (t^ — 1) terms. 

The integrals to compute for the solutions of second order variational 
method are the following 

j if - lfQ^{t)Q^{t)Qk{t)dt j if - lfP,{t)Q^{t)Qu{t)dt 

j - lfP,{t)Pj{t)Qu{t)dt j {e - lfP,{t)Pj{t)Pk{t)dt 

They are all of the form given by Theorem |5j We already know that the first 
one, the third one and the last one satisfy the condition of Theorem |5| so 
they all belong to 

C t, arctanh 1^^^ , In (t^ - l) 
For the second one, we compute the coefficient in a of the residue. This gives 



1 



coeff(Res - lfP,{t)Q,{t)Qu{t)dt,a) = ]-t,tk [ Pi{t)P^{t)Pk{t){e - ifdt 

-1 

which equals to zero because it corresponds to the condition in for the first 
one. So the residue condition is also satisfied, and then thanks to Theorem |5| 
it also belongs to 



C 



t, arctanh ^-^ , In (t^ — l) 
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3.4- Applications 

3. 4-1- The non diagonalizable case 

Let us look first at this diagonalizability condition. We will see it is in 
fact not a strong condition, because if the potential is integrable at order 
one, there is only one possibility for which the Hessian matrix can be non 
diagonalizable, only for the eigenvalue —1. We do not make a complete 
analysis of this case at order 2 because it is not possible to make an efficient 
reduction in this case to produce a nice criterion and this case is very rare 
in practice (still in this case, the Theorem [2] gives integrability constraints 
on a subsystem, but this constraint is not a priori optimal), so it is probably 
more adapted to make an analysis case by case directly in applications. 

Theorem 14. We consider the equation 

X{t) = ^,AX{t) (5) 

with A a matrix in Jordan form. Then this equation has a virtually abelian 
Galois group if and only if 

. Ai^, = (p^-^)ip^+^) yj^th p,eN 

• A is diagonal except maybe for eigenvalue —1 for which the Jordan 
blocks should have a size less than 2. 

Proof. In the non diagonalizable case, the non homogeneous part of the 
variational equation correspond to terms outside the diagonal in the matrix 
A. If the Jordan block have size less tan 2, the equation ^ can be rewritten 

X = -^diX + -^Y Y = -^JiY (6) 

We have in particular that Y satisfy the homogeneous equation for the same 
eigenvalue, because a Jordan block has the same eigenvalue on the diagonal. 
For the case with bigger Jordan block, we would get even stronger conditions 
because then the equation ^ would be a subsystem of equation In our 
case (except for eigenvalue —1), this analysis is not necessary. We compute 
the solutions and we prove that the following function should be in the Picard 
Vessiot field. 

S = J if -l){Q,{t) + t,aPifdt 
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With Theorem [5] and Lemma [13} we know it is enough to study the sequence 



Si = Res {t^ - l){Qi{t) + eiaPi] 



i=oo 



The interesting term is in a, because the coefficient in is always zero. The 
Mgfun package give us a recurrence and then a closed form for this residue 

Oi — 



(2i + l)(z + l) 



which is never zero, except for i = 0. The case i = necessitate the analysis 
of a Jordan block of size 3 (higher Jordan block size would still have this 
system as a subsystem). This gives the equation 

Again, we compute the solution and we use Theorem |5] to prove non com- 
mutativity of the monodromy (in this case, it is very easy because there 
are no parameters). To conclude, we notice that the Galois group is always 
connected because we only take recursively integrations (no algebraic func- 
tions are involved) and because the Galois group of X = ^^rfjX is always 
connected. 



3.4-2. A useful corollary 

Theorem 15. Let V be a meromorphic homogeneous potential of degree — 1 
in dimension n, c a Darboux point ofV with multiplier —1. We pose Xi i = 
1 . . .n the eigenvalues 0/ V^V"(c) with Ai = 2 (the eigenvalue 2 always appear 
in the spectrum). Suppose that V^V^(c) is diagonalizable and 



{A2,...,A„}c{(2fc-l)(A; + l), keB} 
where B G N such that 



ma,x{B) < max(2min(i?) — 1,0) 



Then the variational equation at order 2 near the homothetic orbit associated 
to c has a virtually abelian Galois group. 
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Remark 1. In practice, this corollary says us that if the eigenvalues of the 
Hessian matrix have all an even index and are sufficiently near to each other, 
then the system is always integrable at order 2 without any additional con- 
ditions. Moreover, this theorem is in some sense "optimal", because if it 
is not satisfied, then there will be strong additional integrability constraint 
(of codimension at least 1) for integrability. It allows also to have a strong 
intuition about what will be the easy and the hard cases in proving non 
integrability of a particular problem depending on parameters. 

Proof. We have the Euler relation due to homogeneity and the Darboux 
point condition 



1=1 



By derivating the Euler relation in qj, we get 



With the Darboux point relation, this implies that c is an eigenvector with 
eigenvalue 2. Let us note Xi = c, X2, . . . , Xn a basis of eigenvectors of 
V^V^(c). We will first prove that 

D^{V){c).{Xi,Xa,Xh) = \/a^b 

We derive the Euler relation Q two times and evaluate on c 

^ dqidqjdqk Oqjdqk 
We multiply each line with index j by {Xa)j and we sum over the j 



with Xa the eigenvalue associated to Xa and using the fact that Xa is an 
eigenvector of V^\^(c). We then multiply each line with index k by {Xi,)k 
and we sum over the k 

n n n p. 
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thanks to orthogonality. This is the expression of D^(y){c).{Xi,Xa, Xf,). Let 
us now use Theorem HI We first remark that all invoked index of table A are 
even. Moreover, we have that for three even index i,j,k, if max(i,j, fc) < 
2 min(z, j, k) — 2, then 

A-i^i^i, Ajjj, ^i,j,j'> ^i,k,ki Aj^i^i, AjP;)^^ Af^^i^i, Aj^ j j, Aij^)^ 1 

We also have Ao,o,o = 1- So, if the eigenvalues of V^V^(c) satisfy 

{Ai, A2, . . . , AJ C {(2A; - l)(fc + 1), A; G S} 
where i? C N and such that 



max(i?) < max(2 min(i?) — 1,0) 



then the system is integrable at order 2. Still, knowing that the eigenvalue 
2 always appear in the spectrum of V^V^(c), this would be useless. But 
we know that c is always an eigenvector with eigenvalue 2, and we have 
that all the possible conditions linked to this eigenvector are of the form 
D^(V^)(c).(Xi, Xq, Xfe) = 0. These are automatically satisfied for a ^ h. For 
a = b, we have that v42,i,j = 1 Vz 7^ 1, so the only possible problem would 
be if Xa has the eigenvalue 0, but in this particular case, we also get with 
equation 



n n n 

because A^ = 0. So the possible integrability conditions involving the eigen- 
vector Xi = c are always satisfied, and thus we can remove one time the 
eigenvalue 2 from B. This gives the theorem. 

3.5. Study of the Galois group in the integrable case 

We now prove Theorem |3} analyzing more precisely the Galois group in 
the case where the integrability conditions of Theorem |5] are satisfied. We 
will see that in fact the Galois group almost never grows, and the Galois 



group can be in fact precisely computed thanks to Lemma 13 



Proof. In the integrable case, the variational equations at order 1 involve 
the functions Pi,Qi which are in C(t,arctanh (1)) (after variable changes). 
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The only non multivalued function Q is the function Qo; but the eigenvalue 
2 is always in the spectrum, and so the Galois group is always C. 
At order 2, using Theorem [5} we already know that the solutions are in 
C [t, arctanh (y) ? lii (^^ ~ 1)] and we know a condition for which the term 
in In {t^ — 1) do not appear. Thanks to Lemma 13, we know that this log- 
arithmic term can appear only if Sij^k is a non zero constant (independent 
of a because we suppose that the second order variational equations have a 
virtually abelian Galois group). Let us prove that 

eval(5i J, fc, a = 0) = V i, j, A; G N* (9) 

We only need to use the recurrence Q for Sij^k- To prove that this sequence 
is zero, we then only need to prove it vanishes on the boundary, and here it 
comes down to the case i = j = 1 (because of the invariance by permutation 
and the recurrence on the plane i = 1). Mgfun then prove that 

eval(Si,i,fc, a = 0) = \f keW 

With recurrence (|4]), we get the property 

Let us look now to the case where one of the index is zero. We need to 
study 

Sl^ = Res (t' - l)g,(t)g,(t) ^ j^gg (^2 _ i)tg.(t)g.(t) 

'■' t=oo " t=oo 

We also prove they vanish for a = using recurrence. If two index are zero, 
then we need to study 

Si = Res Qi{t) = Res tQi{t) Sf = Res t'^Qi{t) 

t=oo t=oo t=oo 

All these sequences are zero except for z = 1 for which = — |. Eventu- 
ally, in the case where all index are zero, we need to compute the following 
integrals 

1 ft f f 

-dt / dt / dt / dt 



t'-l J t'-l 

The second and the forth integral have a term in In (t^ — 1). This gives 
Theorem |3l 
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Remark 2. The computation of the sequence eval {Sij^k, a — 0) has in fact 

no sense when Sij^k depend on a. Indeed a corresponds to the multivalu- 
ation Qi. If eval (S'jj^fc, a = 0) = and Sij^k depend on a, then when we 
replace a by a + 1, we get eval {Sij^k, « = 0) 7^ and so this vanishing term 
correspond in fact only to a convention taken for the Q^. Still the convention 
is well chosen, because it allows to study eval a = 0) without making 

a distinction between integrable cases and non integrable cases (in particu- 
lar, eval {Sij^k, = 0) is almost always zero for all values, and this property 
allows a much faster proof than in previous sections). 
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